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Droplet  Sizing  Using  the  Shifrin  Inversion 

R.  ALBERT  and  P.  V.  FARRELL 

Department  of  Mechanical  Engineering 
University  ot  Wisconsin-Madison 
Madison,  Wisconsin  53706 


Abstract 

A  method  for  measuring  droplet  size  distributions  was 
investigated  with  an  emphasis  on  limitations  related  to 
measurements  in  real  spray  environments  The  method  stores  a 
photographic  image  of  a  plane  of  droplets  within  a  spray  and  is 
capable  of  evaluating  particle  size  distributions  within  the  spray, 
one  small  region  at  a  time.  The  method  complement  droplet 
velocity  measurements  made  using  Particle  Image  Velocimetiy.  In 
a  typical  experiment,  a  plane  of  the  spray  was  illuminated  by  a  laser 
light  sheet  and  photographed.  After  processing,  a  :mall  laser  beam 
scanned  the  film  and  a  diffraction  pattern  was  generated  for  each 
region  illuminated  by  the  small  laser  The  diffraction  pattern  was 
inverted  using  a  Shifrin  inversion  to  solve  for  the  particle  size 
distribution  within  the  illuminated  region 

When  applied  to  real  particle  images  for  either 
monodisperse  or  polydisperse  particle  distributions,  the  results  of 
the  Shifrin  inversion  are  affected  by  constraints  on  the  range  and 
quality  of  the  intensity  data  Some  of  these  constraints  are  imposed 
by  hardware  limitations  and  some  are  imposed  by  the  coherent 
reconstruction  system  itself.  In  this  paper  we  will  discuss  some  of 
these  limitations  and  indicate  one  approach  which  seems  to  allow 
for  improved  inversions  using  data  signal  processing. 

Experimental  data  for  use  with  the  measurement  method 
was  generated  using  water  doped  with  2',7’  dichlorofluorescein 
from  an  atomizer 

Nomenclature 

d  =  particle  diameter 

Iq  =  intensity  at  an  angle  0  from  the  optical  axis 

k  =  wave  number.  2rt/>. 

n(x)  =  particle  distribution  function 

x  =  particle  size  parameter,  ndfk 

x,  =  limit  of  integration  for  size  parameter  in  eqn  (1) 

x2  =  limit  of  integration  for  size  parameter  in  eqn  (1) 

y  =  particle  parameter,  x8 

z  =  distance  from  particle  image  to  image  plane  for  interrogation 
X  =  wavelength  of  light 
0  =  angle  from  optical  axis 


Introduction 

A  significant  number  of  two-phase  flow  problems  are  of 
short  duration,  or  exhibit  important  transient  characteristics  upon 
starting  or  stopping,  or  are  in  a  turbulent  flow  field  Any  of  these 
conditions  renders  many  of  the  details  of  the  flow  unsteady 
Measurement  of  the  spatially  resolved  transient  behavior  of  the  flow- 
field  may  provide  better  understanding  of  two-phase  flow  behavior 

Current  experimental  techniques  for  measuring  particle  size 
and  velocity  have  been  carefully  developed  and  appear  to  be  mature 
techniques.  These  include  Fraunhofer  diffraction  measurements 
(Swithenbank,  et  al.,  1976  and  Yule,  et  al.,  1982),  interferometric 
intensity  type  measurements  (Fanner,  1972  and  Bachalo,  1980), 
and  phase  Doppler  type  measurements  (Bachalo  and  Houser, 
1984).  Comparisons  of  these  techniques  are  given  by  Jackson  and 
Samuelson  (1987)  and  Dodge,  et  al  (1987)  and  recent  summaries 
are  given  by  Hirleman,  et  al.  (1990) 

These  techniques  make  measurements  at  a  single  point  or 
measurements  averaged  over  some  small  area  of  the  flow.  They 
can  make  a  series  of  measurements  at  a  given  location  over  time 
For  transient  flows,  particularly  ones  with  significant  turbulence  and 
other  stochastic  characteristics,  a  collection  of  single  point  time 
records  spanning  some  spatial  region  of  interest  is  difficult  to 
interpret.  The  various  spatial  locations  would  be  sampled  at 
different  discrete  times,  for  what  may  amount  to  different  flow- 
conditions  For  many  transient  multi  phase  flow  applications, 
measurements  which  span  a  spatial  area  or  volume  at  a  series  of 
discrete  times  may  be  of  greater  utility  in  developing  an 
understanding  of  the  flow  as  well  as  for  development  and 
refinement  of  flow  models. 

Measurement  Technique  Background 

The  basis  of  the  method  for  measuring  particle  size  and 
velocity  is  similar  to  the  developing  technique  of  particle  image 
velocimetry  (PIV)  A  brief  review  of  the  concepts  necessary  for 
this  measurement  technique  will  be  presented 

Particle  image  velocimetry  has  been  described  in  a  complete 
review  by  Adrian  (1991)  The  general  technique  of  PIV  involves  a 
multiple  exposure  photograph  of  a  flow  containing  particles  or 
droplets  A  photographic  image  of  the  particles  is  obtained  for  a 
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plane  of  particles  which  has  dimensions  of  the  height  of  an 
illuminating  laser  sheet,  the  width  of  the  image  plane,  and  the 
thickness  of  the  laser  sheet  Typical  dimensions  are  3  cm  high  by  5 
cm  wide  by  200  pm  thick.  The  light  source  is  usually  controlled  to 
allow  two  exposures  of  the  particle  field  to  be  recorded  on  the  film. 
The  time  between  the  exposures  is  controlled  and  is  presumed 
known.  Particle  velocity  can  be  obtained  by  measuring  the 
magnitude  and  direction  of  the  displacement  of  a  particle  between 
exposures. 

After  development,  the  film  is  illuminated  with  an 
unexpanded  laser  beam  Sequential  illumination  of  portions  of  the 
film  provides  information  regarding  particle  displacement  in  the 
region  illuminated.  Each  particle  image  generates  a  diffraction 
pattern  in  the  far  field  (far  from  the  film)  The  combination  of  the 
many  diffraction  patterns  from  many  particle  images  illuminated  by 
the  laser  beam  produces  an  interference  pattern,  known  as  Young's 
fringes  The  low  frequency  modulation  of  the  fringe  pattern  is 
known  as  the  diffraction  halo. 

Particle  size  information  is  available  from  the  shape  of  the 
diffraction  halo  which  modulates  the  Young’s  fringe  pattern.  Since 
the  Fourier  transform  or  far-field  diffraction  operator  is  linear,  the 
effects  of  many  particles  may  be  interpreted  through  superposition. 
The  diffraction  halo  is  the  same  diffraction  pattern  for  near  forward 
scattering  from  particles  which  is  tbserved  for  diffraction  based 
particle  sizing  techniques  (Swithenbar  k,  ct  al.,  1976). 

Inversion  of  the  scattered  intensity  data,  minus  the  higher 
frequency  Young’s  fringe  data,  can  provide  a  particle  size 
distribution  of  the  ensemble  of  particle  images  which  appear  within 
the  interrogating  laser  beam  It  may  be  noted  that  in  dilute  regions 
of  a  spray,  the  small  sample  volume  interrogated  may  cover  a  few 
droplets,  not  enough  to  be  reliably  predicted  with  a  continuous  size 
distribution.  Thus  in  inverting  the  diffraction  halo  intensity  pattern 
to  recover  a  particle  size  distribution,  it  is  unlikely  that  an  a  priori 
selected  distribution  function  will  be  successful. 

For  particle  size  measurements  a  TV  camera  is  used  to 
capture  the  diffraction  pattern  and  Young's  fringe  pattern  at  the 
back  focal  plane  of  the  lens,  as  shown  in  Fig.  la,  along  with  a 
typical  photographic  set-up  in  Figs,  lb  and  Ic.  The  use  of  a  TV 
camera  as  a  detector  has  the  advantage  of  simultaneous  use  for  PTV 
applications  and  selectivity  of  geometry  of  data  collection.  With 
the  camera,  information  in  the  intensity  profile  can  be  collected  at 
different  sector  sizes  and  at  different  angles  within  each  sector.  A 
disadvantage  of  using  a  TV  camera  is  the  limited  dynamic  range  of 
the  attached  frame  grabber,  which  in  our  case  is  8  bits  for  each 
pixel  In  our  analysis  some  high  intensity  information  in  the  small 
diffraction  angle  region  was  sacrificed  (overexposed  or  saturated) 
to  pick  up  information  at  higher  scattering  angles. 

The  intensity  pattern  captured  by  the  TV  camera  can  be 
sampled  in  a  variety  of  ways,  using  line  scans  of  the  image  in 
selected  directions,  integrating  rings  at  specified  radii  to  imitate  ring 
detectors  (Kouzalis,  et  al,  1988),  or  integrating  over  selected 
sectors  of  rings.  In  the  experiments  discussed  here,  little  difference 
was  seen  in  the  results  for  line  scans  and  sectors  integrated  over  as 
much  as  45°  Since  line  scans  are  much  faster,  these  are  used  for 
the  results. 
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Figure  la:  Interrogation  system  for  PIV  and  panicle  sizing  method 


(esi  ceil 


cylindrical  lenses 

O 

nozzle 

Laser 

n  . . n 

11  u 

& 

35  mm  camera 


Figure  lb.  Experimental  arrangement  for  spray  photography,  plan 

view 


Figure  lc  Experimental  arrangement  for  spray  photography, 
elevation 

The  diffraction  halo  due  to  panicle  scattering  is  the  low 
frequency  portion  of  the  Young's  fringe  pattern,  which  is  given  by 
the  Airy  formula 


h  )i!20’>*2nU><£c  (1) 

k  V  <  JC, 

where  0  is  an  angle  measured  from  the  center  of  the  diffraction 
pattern,  Ifl  is  the  normalized  scattered  intensity  at  an  angle  0,  k  = 
2 tc/X,  A.  is  the  wavelength  of  light,  z  is  the  distance  from  the 
diffracting  object  to  the  viewing  plane,  x  is  the  particle  size 
parameter  x  =  nd/k,  Xj  and  x2  are  limits  on  the  expected  particle 
size  parameter,  n(x)  is  the  particle  size  distribution,  0  »  sin  0,  y  = 
x9,  and  J)  is  a  Bessel  function  of  the  first  land  of  order  unity  For 
some  distributions  n(x)  may  be  assumed  to  vanish  at  x.  and  x2, 
giving 
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Equation  (2)  may  be  considered  a  first-kind  Fredholm 
integral  equation  An  inversion  producing  n(x),  known  as  a  Shifrin 
inversion,  is  possible  (Fymat  and  Mease,  1978)  resulting  in 
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where  Yj  is  a  Bessel  function  of  the  second  kind  of  order  1 

This  type  of  data  inversion  allows  flexibility  in  the  type  of 
particl*  size  distribution  acquired,  si"'-"  no  distribution  is  assume:1 
A  major  disadvantage  of  this  integral  inversion  scheme  is  that  the 
measured  data,  I  ,  must  be  differentiated.  Alternative  inversion 
methods  using  the  Shifrin  inversion  in  other  forms  (Fymat  and 
Mease,  1978  and  Coston  and  George,  1991)  and  the  use  of  other 
inversion  schemes  (Hirleman,  1991)  have  been  suggested.  Some  of 
these  inversions  are  aimed  primarily  at  ring  detector  geometry,  but 
others  which  may  prove  useful  in  this  application  have  not  yet  been 
implemented  It  should  be  noted  that  the  result  of  the  calculation 
from  Eqn  3  is  a  relative  particle  size  distribution,  and  actual 
particle  flux  or  mass  flux  of  particles  is  not  a  result  of  Eqn  3. 
Estimates  of  total  particle  flux  for  diffraction  based  methods  are 
usually  made  by  measurement  of  the  attenuation  of  the 
interrogating  laser  beam  relative  to  a  correlation  for  extinction  as  a 
function  of  flux 

The  technique  outlined  has  been  tested  on  modeled  and 
experimental  single  apertures,  a  variety  of  modeled  polydisperse 
distributions,  a  fixed  calibration  reticule,  and  a  monodisperse 
droplet  stream  (Farrell,  1991).  The  results  of  this  previous  work 
indicated  that  the  inversion  scheme  worked  well  for  both 
monodisperse  and  polydisperse  size  distributions  using  modeled 
intensity  data  Using  experimental  data  from  single  apertures,  the 
method  was  able  to  accurately  measure  monodisperse  distributions. 
Using  a  reticule  with  a  Rosin-Rammler  distribution  of  plated  spots, 
the  method  was  able  to  accurately  measure  the  distribution.  An 
example  of  the  size  distribution  measured  for  the  Rosin-Rammler 
distribution  of  the  spots  on  a  reticule  is  shown  in  Fig  2.  The  results 
for  the  distribution  agree  fairly  well  with  the  "given"  distribution  for 
the  images  In  the  previous  work,  the  method  had  not  been 
extensively  tested  in  spray  fields  subject  to  laser  speckle  and  other 
noise  sources 

When  the  experimental  method  described  was  applied  to  a 
spray,  several  problems  were  observed  which  were  not  clearly 
eviden*  in  the  single  aperture  and  polydisperse  reticule  evaluations 
For  a  typical  application,  a  large  spray  field  is  photographed,  either 
with  single  or  double  exposure,  and  interrogated  on  a  scale  of  about 
0  5  mm  For  many  regions  of  a  spray,  this  geometric  scale  will 
include  a  small  number  of  droplets  In  these  regions,  system 
performance  for  monodisperse  or  multi  mode  size  distributions  is 
important. 
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Figure  2  Experimental  distribution  profile  for  Rosm-Rammier 
distribution,  compared  with  data  supplied  with  reticule 


Particle  Sizing  Method  Improvements 

The  panicle  sizing  method  may  be  charactenzed  using 
calculated  intensity  distributions  for  monodisperse  or  polydisperse 
panicle  distributions  Experimental  monodisperse  and  polydisperse 
distributions  can  also  be  evaluated,  although  quantitative 
comparison  for  polydisperse  distributions  is  difficult  since  the  actual 
distribution  is  rarely  known 

A  calculated  intensity  profile  for  a  50  pm  panicle  is  shown 
in  Fig  3.  This  figure  shows  a  slice  through  the  cross  section  of  the 
Airy  pattern  described  in  Eqn.  (2)  The  numerical  intensity  pattern 
has  been  truncated  in  the  figure  to  expand  the  intensity  scale  at  the 
lower  intensity  values,  but  the  complete  intensity  profile  is  used  in 
the  inversion  calculations  Using  the  Shifrin  inversion,  the  area 
weighted  number  distribution,  x2n(x),  calculated  using  the  Shifrin 
inversion  is  shown  in  Fig  4  A  50  pm  aperture  intensity  profile  was 
measured  experimentally,  and  the  area  weighted  size  distribution 
was  calculated.  These  profiles  are  also  shown  in  Fig  3  and  Fig  4 
Note  that  the  experimental  intensity  profile  in  Fig  3  has  saturated 
the  detector  (TV  camera)  for  the  first  two  peaks  As  indicated  in 
Fig  4,  the  actual  aperture  diameter  was  slightly  larger  than  50  pm 
It  was  measured  to  be  about  53  pm. 

In  comparison  to  the  numerical  intensity  distribution  in  Figs 
3  and  4,  the  experimental  distribution  shows  greater  noise  and  a 
limited  range  of  6  over  which  the  intensity  signal  is  above  the  noise 
floor.  As  a  result,  the  integration  in  the  Shifrin  inversion  is  carried 
out  over  a  smaller  range  of  significant  6  values,  and  the  resulting 
distribution  in  Fig.  4  shows  broader  peaks  and  some  "ringing"  near 
the  50  pm  peak  For  this  particular  experiment,  the  noise  in  the 
intensity  profile  is  primarily  background  optical  noise  (e  g  laser 
speckle)  so  camera  and  digitizer  noise  are  not  the  major 
contributors.  If  the  ripple  in  the  distribution  in  Fig  4  is  considered 
to  be  random  noise,  an  RMS  noise  value  can  be  calculated  for  each 
of  the  profiles  in  Fig  4  The  numerical  profile  has  an  RMS  noise 
figure  of  0  761%  of  the  peak  signal,  while  the  experimental  profile 
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has  an  RMS  noise  figure  of  3  465%  of  the  peak  signal  Clearly  the 
noise  in  these  distributions  is  not  random,  and  the  RMS  noise  does 
not  accurately  characterize  the  ability  to  identify  distribution  peaks, 
but  this  type  of  calculation  gives  a  numerical  comparison  of  the 
ripple  of  the  distribution  function 
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Figure  3:  Intensity  profiles  for  numerical  and  experimental  data 


figs  5  and  6  in  which  numerically  generated  intensity  profiles  for  a 
50  pm  panicle  have  been  truncated  to  smaller  total  sampling  angles 
In  Fig  5,  the  number  of  sample  values  is  changed  while  the  A0  is 
held  constant,  and  in  Fig  6,  the  number  of  sample  values  is  held 
constant  and  A0  is  allowed  to  vary  In  either  case,  the  overall  angle 
for  sampling,  0^,  has  a  significant  effect  on  the  size  distribution 
pattern  The  RMS  noise  in  the  size  distribution,  calculated  as  for 
Fig  4,  ranges  from  76%  to  11  44%  in  Fig.  5  and  from  76%  to 
11.29%  in  Fig  6. 

When  using  intensity  profiles  for  monodisperse 
distributions,  the  required  range  of  0  is  significantly  larger  than  for 
the  polydisperse  case  This  increase  in  required  effective  sampling 
angle  for  monodisperse  distributions  has  been  previously  noted  by 
Fymat  and  Mease  (1978)  In  the  current  measurement  method,  a 
large  portion  of  the  particle  image  photograph  lies  in  a  region 
where  a  small  number  of  particles  will  be  sampled  by  the 
unexpanded  laser  beam  at  any  one  time  Thus,  performance  on 
monodisperse  distributions  or  small  numbers  of  panicle  sizes  is 
important  in  the  overall  performance  of  the  system 


Figure  5;  Numerical  data  with  variation  of  0mav  A0  constant 


Figure  4:  Comparison  of  area  weighted  size  distributions  for 
numerical  and  experimental  data 

Reducing  the  range  of  0  over  which  the  numerical  intensity 
profile  is  sampled  or  has  significant  content  produces  characteristics 
similar  to  the  experimental  aperture  results,  indicating  that  the 
limitation  to  relatively  small  0  values  has  a  significant  impact  on  the 
calculated  size  distribution  and  the  distribution  pattern.  The  effect 
of  reduced  total  angle  over  which  data  can  be  collected  is  shown  in 
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Figure  6:  Numerical  data  with  variation  of  0max,  Ne  constant 

Several  approaches  to  "extrapolate"  the  measured  data  to 
larger  0  values  for  monodisperse  distributions  have  been  tried,  and 
one  shows  some  promise  Theoretical  particle  size  distributions 
produce  intensity  distributions  which  are  essentially  integrals  of 
Bessel  functions  (Eqn.  1).  A  part  of  the  kerne,  of  the  integral  used 
in  the  Shifrin  inversion  is  the  ter... 
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an  effect,  giving  larger  peak  values  for  the  area  weighted  vze 
distribution,  and  less  npple  in  the  size  distribution  An  example  of 
the  results  of  this  type  of  data  extrapolation  is  shown  in  F:g  7 
This  figure  shows  the  results  of  two  inversions  from 
experimentally  derived  intensity  distribution  for  a  nomine. y  50  um 
aperture  One  of  the  inversions,  labeled  with  the  stars,  is  the 
inversion  from  the  original  data  set,  without  extrapolation  o; 
filrering.  The  second  profile,  labeled  with  the  circles,  is  from  tne 
transformed,  interpolated,  filtered,  and  inverse  transformer!  data 
The  two  profiles  are  clearly  similar,  but  the  extrapolated  one 
exhibits  a  slightly  wider  profile  at  the  50  ur  particle  -...e  .  c 
significantly  reduced  npple  elsewhere  Using  the  RMS  noise  a«  an 
estimate  of  the  ripple  in  the  size  distribution,  the  "original  data"  in 
Fig  7  has  an  RMS  noise  of  3  46°  i  ar.d  the  "extrapolated  data"  has 
an  RMS  noise  of  1  82% 
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which  for  numerically  generated  monodisperse  size  distributions  is 
essentially  a  cosine  wave  In  order  to  extrapolate  the  range  of  0 
over  which  this  term  is  available  (non-zero  and  not  dominated  by 
noisej,  the  existing  portion  of  the  kernel  is  treated  in  a  manner 
similar  to  those  used  for  calculations  of  super  resolution  in  optics. 
The  kernel  cited  in  Eqn  4,  as  a  function  of  0,  is  Fourier 
transformed  using  an  FFT  routine  In  the  Fourier  plane,  the 
available  frequency  values  may  be  used  to  generate  more  values  by 
interpolating  between  calculated  frequency  values  In  our 
implementation,  using  a  trapezoidal  rule,  this  technique  effectively 
doubles  the  number  of  frequency  samples.  Upon  inverse  Fourier 
transformation,  the  reduction  of  Af  by  a  factor  of  2  leads  to  an 
increase  in  the  total  9  range  of  the  function  by  a  factor  of  2  This 
method  assumes  that  all  significant  frequency  content  of  the  kernel 
is  already  present  in  the  original  data  The  result  represents 
extrapolated  data  which  is  consistent  with  the  original  values  in 
terms  of  frequency  content  In  an  effort  to  filter  the  high  frequency 
noise  typically  present  on  experimental  signals,  we  also  window  the 
frequency  spectrum  before  inverse  transformation  using  a 
Blackman-Harris  type  window  After  this  extrapolation  procedure, 
the  newly  lengthened  kernel  is  used  in  the  Shifrin  inversion  to 
determine  the  particle  size  distribution 

This  technique  has  almost  no  effect  on  numerically  derived 
intensity  distributions,  since  typically  any  range  of  0  values  may  be 
selected  when  the  intensity  distribution  is  generated  For 
monodisperse  experimental  data,  the  extrapolation  of  the  data  has 
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Figure  7  Comparison  of  original  experimental  data  with 
extrapolated  data 

Spray  System 

Experimental  results  from  a  water  spray  are  shown  to 
indicate  how  the  method  performs  on  an  actual  spray  Since  the 
size  distribution  in  the  regions  interrogated  by  the  laser  is  not 
known,  only  qualitative  evaluations  of  the  system  performance  is 
possible  under  these  conditions. 

A  spray  chamber  was  designed  and  bui't  of  acrylic  with 
three  windows  to  serve  as  entry  and  exit  points  for  the  laser  sheet 
and  as  a  window  for  photography  A  Spraying  Systems  nozzle, 
which  was  an  external  mix,  round  spray,  full  cone  atomizer  (I  4J. 
64SS  air  cap,  1650SS  fluid  cap.  with  No  12810  liquid  shut  off 
needle)  was  used  The  chamber  allowed  both  horizontal  and 
vertical  nozzle  translation  relative  to  the  laser  sheet  location  The 
air  was  supplied  under  pressure  from  the  house  supply  directly  to 
the  nozzle  and  the  liquid  was  siphon  fed  from  a  500  ml  beaker 

A  doubled  Nd  Yag  laser  supplied  illumination  of  the  sprav 
by  providing  high  intensity  light  in  a  pulse  duration  of  about  6  ns 
operating  at  a  wavelength  of  532  nm  Two  cylindrical  lertses 
shaped  the  laser  beam  The  light  was  expanded  by  the  first  lens  to  a 
height  of  40  mm  at  the  test  chamber  The  second  cylindrical  lens 
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decreased  the  7  mm  width  of  the  beam  to  about  1  mm  at  the 
location  of  the  spray  A  sketch  of  the  system  is  shown  in  Figs,  lb 
and  lc 

A  35  mm  camera  with  a  105  mm  focal  length  lens  with  a 
focusing  tube  was  used  to  take  the  photographs  of  the  spray  A 
magnification  of  11  was  chosen  for  simplicity  The  imaged 
droplets  at  90°  scatter  appeared  as  a  front  and  a  back  arc,  or  in 
some  cases  as  a  disk  For  the  analysis  scheme,  filled-in  images  were 
needed.  To  solve  this  problem  2', 7'  Dichlorofluorescein  was  used, 
to  dope  the  water  spray  to  provide  fluorescence  from  the  532  run 
source  The  sprayed  solution  contained  1  g  of  the  chemical  for 
every  200  ml  of  water  and  the  remaining  solids  were  filtered  with 
Whatman  #1  filter  paper  A  hydrometer  determined  the  specific 
gravity  of  the  solution  to  be  the  same  as  that  of  water  The 
Dichlorofluorescein  is  listed  as  a  low  hazard  material  with  usual 
industrial  handling.  A  wet  vacuum,  equipped  with  a  foam  filter 
sleeve,  was  connected  to  the  bottom  of  the  spray  chamber  and  the 
exhaust  was  vented  directly  into  the  outgoing  building  exhaust 

The  parameters  which  could  be  varied  during  the 
experiment  were  the  supply  air  pressure  (which  correspondingly 
altered  the  liquid  flow  rate)  and  the  horizontal  and  vertical  nozzle 
positions  For  this  paper,  one  setting  of  these  values  was  used, 
since  the  emphasis  was  on  the  inversion  scheme 

In  the  region  of  interrogation,  the  droplet  images  on  the 
film  typically  range  from  15  to  70  jim.  and  are  relatively  well 
dispersed  Original  and  extrapolated  data  results  are  plotted  in 
Figs  8  and  9  for  samples  taken  at  two  different  locations  of  the 
water  spray  from  the  air-assisted  atomizer  These  locations  were 
chosen  arbitrarily  to  indicate  the  general  types  of  particle 
distributions  available,  and  are  not  intended  to  comprehensively 
describe  the  spray  itself  In  each  figure,  a  plot  of  the  size 
distribution  from  unextrapolated  data  is  shown  in  comparison  with 
the  distribution  using  extrapolated  and  filtered  data  Note  that  for 
these  figures  the  particle  size  distribution  is  plotted,  not  the  area 
weighted  distribution  The  regions  sampled  have  a  few'  to  several 
particles  within  the  region,  and  the  size  distributions  should  reflect 
that  kind  of  multi  mode  character  The  extrapolated  data  appear  to 
have  somewhat  reduced  oscillation  in  the  n(x)  profiles  with  respect 
to  the  unextrapolated  versions  Each  sampled  region  has  only  a  few 
particles  in  it,  so  each  peak  in  the  size  distribution  may  be 
interpreted  as  one  or  more  panicles  at  that  size 

The  size  distributions  shown  in  Figs  8  and  9  are  not  as  well 
defined  as  those  for  the  single  apertures,  in  Fig  7,  for  example 
Figure  8  shows  a  large  peak  at  about  20  p.m  and  smaller  peaks  at 
about  30  pm  and  70  pm,  all  of  which  are  likely  to  represent  single 
particles  of  these  respective  sizes  The  remainder  of  the  distribution 
ir.  Fig  8  shows  remnants  of  the  ripple  seen  on  other  e;;amp!e„.  but 
looking  at  the  "extrapolated"  results,  the  remaining  peaks  do  not 
appear  to  rise  above  the  background  ripple  Using  the  RMS  noise 
calculation  used  in  Fig  7  comparisons,  the  extrapolated  size 
distribution  for  Fig  8  exhibits  an  RMS  noise  value  of  3  95%  of  the 
peak  value,  while  the  original  data  has  an  RMS  noise  figure  of 
5  66%  of  the  peak  value  In  Fig  9,  large  values  of  the  distribution 
function  are  evident  at  about  15  pm  and  again  at  about  25  pm,  39 
pm,  and  55  pm  Each  of  these  peaks  likely  represents  a  single 
droplet  of  that  size  The  remainder  of  the  distribution,  at  larger 
droplet  sizes,  appears  to  be  primarily  background  ripple  The  RMS 
noise  figure  for  the  extrapolated  data  in  Fig  9  is  1 1  5%  of  the  peak 
value,  while  the  original  data  has  a  noise  figure  of  16  5%  of  the 
peak  value  In  both  of  these  cases,  the  use  of  the  extrapolation 


scheme  appears  to  reduce  the  magnitude  of  the  background  npple 
due  to  limited  angular  data  range,  making  the  distribution  peaks 
more  readily  identifiable 
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Figure  8  Sample  size  distribution  from  spray 


diameter  (urn) 

Figure  9  Sample  size  distribution  from  spray 
Conclusions 

This  paper  discusses  a  technique  for  droplet  sizing  in  a 
polydisperse  atomizer  spray  The  major  focus  of  the  paper  is  to 
indicate  how  the  current  inversion  scheme's  level  of  performance  is 
affected  by  limited  data  collection  range  caused  by  experimental 
noise  A  method  is  proposed  for  extrapolation  of  the  data  to  larger 
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sampling  angles  through  interpolation  in  the  frequency  domain 
Some  examples  of  the  performance  of  the  system  are  presented  for 
single  apertures  and  for  portions  of  a  spray  On  single  aperture 
data,  the  "noise"  in  the  size  distribution  interpreted  as  random  RMS 
noise,  decreases  by  a  factor  of  almost  2,  compared  to  a  similar 
"noise"  measurement  on  unextrap  Mated  data  On  experimental  data 
from  a  water  spray,  size  distributions  are  qualitatively  improved  and 
size  peaks  made  more  identifiable  The  reduction  in  noise,  using 
the  calculated  RMS  noise  as  a  figure  of  merit,  is  smaller  for  the 
spray  data  than  for  the  apertures,  with  typical  reductions  of  about 
1/3  ' 
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